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The in-plane infrared response of the high-T c cuprate superconductors was studied using the 
spin-fermion model, where charged quasiparticles of the copper-oxygen planes are coupled to spin 
fluctuations. First, we analyzed structures of the superconducting-state conductivity reflecting 
the coupling of the quasiparticles to the resonance mode observed by neutron scattering. The 
conductivity a computed with the input spin susceptibility in the simple form of the mode exhibits 
two prominent features: an onset of the real part of a starting around the frequency ujq of the 
mode and a maximum of a related function W(ui), roughly proportional to the second derivative of 
the scattering rate [l/r](w), centered approximately at uj = uia + Ao/h, where Ao is the maximum 
value of the superconducting gap. The two structures are well known from earlier studies. Their 
physical meaning, however, has not been sufficiently elucidated thus far. Our analysis involving 
quasiparticle spectral functions provides a clear interpretation. Second, we explored the role played 
by the spin-fluctuation continuum, whose spectral weight is known to be much larger than the one 
of the mode. We have shown that the experimental spectra of 1/r can be approximately reproduced 
by augmenting the resonant-mode component of the spin susceptibility by a suitable continuum 
component with a considerably higher spectral weight and with a characteristic width of several 
hundreds meV. The computed spectra of f/r display a new structure in the mid- infrared which 
is related to the finite width of the occupied part of the conduction band. Third, we investigated 
the temperature dependence (TD) of a assuming that the normal state spin susceptibility consists 
of an overdamped low energy mode and the continuum component. The differences between the 
experimental normal-state spectra and those of the superconducting state, including some interesting 
effects at higher frequencies, are reasonably well reproduced. Motivated by recent experimental 
(ellipsometric) works by Molegraaf and coworkers [H. J. A. Molegraaf et ai, Science 295, 2239 
(2002).] and Boris and coworkers [A. V. Boris et al, Science 304, 708 (2004).], we further studied 
the TDs of the effective kinetic energy K.E. and of the intraband spectral weight To- Calculations for 
the trivial case of noninteracting quasiparticles in the normal state and a BCS-like superconducting 
state reveal a strong sensisitivity of the TD of Jo to details of the dispersion relation. The TDs 
of K.E. and Io in the interacting case, for the set of the values of the input parameters used 
throughout this work, are similar to those of the trivial case. The physics beyond the changes 
occuring when going from the normal to the superconducting state, however, is shown to be more 
complex, involving, besides the formation of the gap, also a feedback effect of the spin fluctuations 
on the quasiparticles and a significant shift of the chemical potential. 



I. INTRODUCTION 



Models of the electronic structure of the high-T c 
cuprate superconductors (HTCS), where charged quasi- 
particles of the copper-oxygen planes are coupled to spin 
fluctuations (SF) (for review see Refs. 1,2) possess two 
appealing features, (a) They provide a straightforward 
explanation of the symmetry of the order parameter and 
(b) the same value of the coupling constant leads both to 
a fairly good agreement between theory and experiment 
for the normal state and to values of T c of 50 — 100 K. 

On the other hand, these models suffer two important 
deficiencies, (a) They cannot be rigorously derived start- 
ing from a well established microscopic hamiltonian such 
as the hamiltonian of the 2D one-band Hubbard model. 
They can only be motivated by some perturbation expan- 
sions of the latter. Within the so called conserving fluc- 
tuation exchange (FLEX) approximation 3,4 ' 5 ' 6,7 ' 8 ' 9 , e.g., 
an important term in the expression for the quasiparticle 



selfenergy E(k, iw„) reads 
1 3U 2 

gN 2 ^2 XRPA(q,i^m)G(k-q,kj„-ift m ), (1) 

q,iO m 

where iuj n (iO TO ) are the fermion (boson) Matsubara fre- 
quencies, [3 = 1/(/sbT), N is the number of momentum 
points, and U is the Hubbard parameter. Further, G 
is the dressed quasiparticle propagator and Xrpa the 
random-phase-approximation(RPA)-bascd expression for 
the spin susceptibility The expression for the selfenergy 
used within the SF-based models has precisely the same 
structure: 

1 3g 2 

J^— XSF(q,i^m)G(k - q,iw„ - ifi ro ) , (2) 

where g and xsf are the model coupling constant and 
spin susceptibility, respectively. This can be considered 
as a kind of formal justification of the SF-based ap- 
proaches. A typical value of g of 0.5 eV, however, is 
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an order of magnitude lower than that of U resulting 
from first-principles calculations of 5 — 10 eV 10,11 . (b) 
The second problem of the SF-based models concerns the 
choice of the function xsf- In the spirit of the models, 
Xsf should match the true spin susceptibility x acces- 
sible to experimental investigations. Unfortunately, the 
available experimental data are not yet sufficient to de- 
termine x hi the relevant spectral range. This applies 
especially to the normal state (NS) of optimum doped 
and ovcrdoped materials, where the spin fluctuations are 
relatively weak. The form of xsf has thus to be guessed, 
at least to some extent. In the classical series of papers, 
Monthoux et al. 12 ' 13 ' 14 have adopted the Millis-Monien- 
Pincs formula 15 : 



XMMp(q,w) 



XQ 



i + (q-Q) 2 C 2 -W^ s 



(3) 



where Q = {it /a, ir /a) is the antifcrromagnctic wave vec- 
tor, a is the in-plane lattice parameter, £ the antifcr- 
romagnctic coherence length, and u> s f is the parameter 
that specifies the frequency of the relaxation mode (at a 
fixed q, x^j MP (q, w) exhibits a broad maximum centered 
at u> = cj s f[l + (q — Q) 2 £ 2 ], i.e., a relaxation mode). 

The situation is simpler for the superconducting state 
(SCS), where a sharp resonance in x ("magnetic mode") 
centered at q = Q, Q = (ir/a, 7r/a), and lu = ujq, 
luq ~ 40 meV, has been detected in neutron scattering ex- 
periments (see Rcfs. 16,17, 18,19 and references therein). 
For some compounds the shape of the resonance is well 
known and it is thus possible to investigate the coupling 
between the mode and the charge carriers without any ad 
hoc assumptions regarding the form of xsf- Let us em- 
phasize that there must be a coupling between the well 
defined mode and the well defined Bogoliubov quasipar- 
ticles and it should be considered in any theory of super- 
conductivity in the HTCS, in particular, when describing 
single-particle excitations (observed by photoemission) 
and two-particle excitations (observed, among others, by 
optics). Such considerations may help to clarify the role 
played by the SF in the mechanism of superconductiv- 
ity. The single-particle case has been reviewed in detail 
in Refs. 1 and 20. An important observation is that the 
peak-dip-hump structure observed in superconducting- 
state angle-resolved-photoemission spectra (for review 
see Ref. 21) at and around (tv/cl, 0) can be understood 
and described in terms of the above mentioned coupling. 
This result, together with several other findings, suggests 
that the SF are an important player in the mechanism. 
In this work we focus on the influence of the coupling on 
the in-plane infrared response of the HTCS in the super- 
conducting state. 

The SF-based approach has been pioneered in this 
context by Schachingcr, Carbotte, and Marsiglio 22 and 
by Quinlan, Hirschfeld, and Scalapino 23 . The former 
group used the boson spectral density proportional to 
XmmpCi = Q'^- 1 )' m the work of the latter xrpa has 
been used. The real part of the in-plane infrared con- 
ductivity, (Ti(w), resulting from these computations ex- 



hibits the following trends: (i) Below T c and for frequen- 
cies lower than ~ 4A , where A is the maximum value 
of the superconducting gap, it decreases with decreas- 
ing temperature, (ii) It does not possess a true energy 
gap, even at very low temperatures and in the absence of 
any impurity scattering. Instead, at low temperatures, it 
increases gradually with increasing to. (iii) The increase 
becomes steeper above the characteristic spin-fluctuation 
frequency of u) s f = 30mcV 22 or Ao - 2Ao 23 . (iv) A very 
broad maximum appears in the low-temperature spectra 
around ~ 4A . These findings are roughly consistent 
with experimental data 24 ' 25 ' 26 ' 27 ' 28 ' 29 , except for two fea- 
tures. First, the data contain additional sharp structures 
at low frequencies. These are probably related to charge 
inhomogeneities and will not be discussed here. Second, 
the increase of o\ (uj) in the data appears to be less grad- 
ual than in the theoretical spectra. The data rather seem 
to exhibit an onset starting around 300 — 400 cm" 1 . This 
may be related to the presence of a sharp feature in the 
spin-fluctuation spectra below T c , i.e., to the magnetic 
mode, which has not been considered in Refs. 22 and 23. 

The mode has been taken into account by two of the 
authors and M. Cardona (MBC) 30 , who employed a form 
of Xsf reflecting the results of the neutron scattering ex- 
periments: 



XRM(q,^) = 



l 



i + (q - Q) 2 e ujZ-lj 2 - iTcj ' 



(4) 



where T is the broadening of the resonance mode and F 
expresses its "oscillator strength". This study provided 
a tentative interpretation of the onset feature: the final 
states consist, in the first approximation, of two Bogoli- 
ubov quasiparticlcs and the magnetic mode; the mini- 
mum excitation energy is thus Tiloq ps 40 meV ps 320 cm -1 
and the broad maximum around 1000 cm -1 corresponds 
to excitations involving quasiparticles around the saddle 
point and the mode. Carbotte, Schachinger, and Basov 
(CSB) 31 also interpreted the data in terms of the mag- 
netic mode but an inverse strategy has been used. Their 
approach is based on a relation between the conductiv- 
ity and the electron-phonon spectral density a 2 F(uj) that 
applies to the normal state of a weakly coupled isotropic 
electron-phonon system 32 : 



a 2 F(uj) ps W{u) 



£ o^pi d 2 
2tt duj 2 



wRe 



1 



(5) 



where cu p \ is the plasma frequency of the free charge car- 
riers. CSB assumed (and to some extent verified by com- 
puting o~(u>) within Eliashbcrg formalism) that a slightly 
modified equation, 



,) PS W{w) 



2ir duj 2 



loRc 



cr(w) 



(6) 

can be used to obtain the "effective electron-spin- 
fluctuation spectral function" I 2 Xef{^>) of a d-wave 
superconductor 22 . The main point of CSB is that Xe/(w), 
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as determined from the optical data, contains a sharp 
structure that is fairly similar to the one of x"(Q,cj), 
i.e., to the neutron resonance. Based on this finding, 
CSB concluded that the infrared data reflect a coupling 
of the charge carriers to the magnetic mode. The cou- 
pling strength inferred from experiment is found to be 
sufficient to account for the high value of T c . The re- 
lation between <Ji{uj) and the magnetic mode has been 
further explored by Abanov, Chubukov, and Schmalian 
(ACS) 33 ' 34 . For example, they predicted a sharp onset of 
<j\(lu) located at htu = 2A + hu> (see Fig. 2 of Ref. 33). 
The three approaches (MBC, CSB, ACS) differ in their 
emphasis and computational formalism. MBC showed 
that the superconducting-state infrared data can be un- 
derstood in terms of the SF-based model with the spin 
susceptibility in the form of the neutron peak and fo- 
cused on the onset in the conductivity spectra starting, 
according to the model, around too. CSB found an ap- 
proximate solution of the inverse problem, the spectral 
function W(uj), and attributed its main maximum to the 
magnetic mode shifted by Ao towards higher frequencies. 
Finally, ACS concentrated on the contribution to the con- 
ductivity of the so called hot spots, i.e., intercepts of the 
Fermi surface and the antiferromagnetic-Brillouin-zone 
(BZ) boundary, and they argued that this contribution 
is the dominant one. Regarding the formalism, all the 
three groups expressed the quasiparticle selfenergy us- 
ing the SCS version of Eq. 2 and the conductivity using 
the textbook formulas 35 , where the vertex corrections are 
neglected. The differences consist in treating the propa- 
gator G and the susceptibility xsf- In the work of MCB, 
G is replaced with the bare Nambu matrix containing 
an estimated superconducting gap of d x i_ y i symmetry, 
i.e., with the BCS expression. The approach is thus not 
self-consistent. On the other hand, it has two important 
advantages, not shared by the other two approaches, be- 
sides its formal simplicity: (a) It is not restricted to mod- 
els, where the SF are the only cause of superconductivity, 
(b) The full g-dependence of \ is taken into account. The 
same approach has been used in the studies of the single- 
particle properties by Eschrig and Norman and coworkers 
(see Ref. 20 and references therein) and, in the context 
of electron-phonon interaction, by Sandvik, Scalapino, 
and Bickers 36 . CSB obtain the propagator G by solving 
Eliashberg equations with a separable interaction, i.e., in 
a selfconsistcnt way. ACS use an approximate form of 
these equations valid in a region around the hot spots. 
Not only the propagator but also the spin susceptibility 
is treated in a self-consistent way. 

The structures in the conductivity spectra can be rela- 
tively easily related to and understood in terms of those 
of the quasiparticle spectral functions (in the absence of 
vertex corrections). Such an analysis, however, has not 
yet been performed, except for some observations regard- 
ing the contributions of the hot spots 30,34 and the optical 
scattering rate 37 . One of the aims of the present paper is 
to provide a detailed interpretation of the onset starting 
around loq discussed by MBC and of the main maximum 



of the function W(uj) emphasized by CSB in terms of 
the quasiparticle spectral function. Among others, we 
explore the role played by the hot spots. 

The structure of the paper is the following. In Sec. II 
we summarize the basic equations of our approach and 
in Sec. Ill we discuss the input parameters and present 
some computational details. Section IV contains our re- 
sults. We start with the quasiparticle spectral functions, 
that are compared with those computed by Eschrig and 
Norman 20 (Subsec. IV A), and contributions of the indi- 
vidual k-points to the real part of the conductivity (Sub- 
sec. IV B). Subsecs. IV C and IV D deal with the spec- 
tra of the total conductivity, in particular with its onset 
around ujo, and with the structures of the scattering rate 
and the function W(cu), respectively. In Subsec. IVE we 
comment on the role of the SF continuum. Some aspects 
of the temperature dependence (TD) of the spectra, as 
the role of the TD of xsf, and the TD of the integrated 
spectral weight, that has been recently related to changes 
of the effective in- plane kinetic energy 29 ' 38 , are addressed 
in Subsec. IV F. A summary and our conclusions are pre- 
sented in Sec. V. Readers who are interested only in the 
main findings of the paper can skip Sees. II and III and 
some more technical parts of Sec. IV. 

II. THEORETICAL FRAMEWORK 

Our starting point is the Hamiltonian of the spin- 
fermion model: 

H = H + H int , (7) 

where 

#o = X/ kC k« Cka . ( 8 ) 

kev 

#int=.g]Ts(q)-S(-q). (9) 
q 

Here £k is the quasiparticle dispersion, a is the spin index, 
g is the spin-fermion coupling constant, 

s (q) = -7=5 c U^p c ^ ( 10 ) 

is the Fourier component of the electronic spin, r = 
( T ii T 2i T 3) is the vector of the Pauli matrices, and S(q) 
the spin-fluctuation operator. The latter quantity is de- 
fined so that its retarded Green's function is equal to the 
negatively taken spin susceptibility, 

oo 

Xy(q,w) = ~ J Xii(q,*)e ih "dt (11) 

— OO 

with 

X«(q,t) = i0(t) ([«i(q,t),*i(-q,0)]) . (12) 
It is assumed that Xijfau) is isotropic, i.e., Xy (q, w) = 
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A. Quasiparticle selfenergy and spectral function 



with 



We follow here the approach of Ref. 30 and express 
the quasiparticle selfenergy E(k, iE n ) using second order 
perturbation theory (H int being the perturbation) and 
starting from a BCS state with an estimated supercon- 
ducting gap of d x 2_ y 2 symmetry. The resulting formula 
for S(k, iE n ) (2 by 2 matrix) represents a SCS version of 
Eq. 2: 

E(k, iE n ) = 



„ ,, .„ , iE n T + (e k - A*)t3 + A k ri 

Go(k '^ ) = W^^AT- (14) 

Here Go is the bare Nambu Green's function \E n are the 
fermion Matsubara energies (iE n = \huj n ), fx is the chem- 
ical potential, and A k the superconducting gap. The full 
Nambu propagator G is given by the Dyson equation, 
G -1 = Gq 1 — E _1 , and the matrix spectral function 
A(k,E) = -2Im[G ret (k,£)] by 



_l_3^ 

fiN 4 



X] XSF(q, «^m)G (k - q, iE„ 



A(k, E) = -2 Im [G J (k, \E n ) - E(k, LE„)] iE ^ E+iS . 

iWlm) (15) 
Equation 13 can be further rewritten into a form which 
(13) is more suitable for practical computations: 



(16) 



where density. The two terms in the bracket on the right hand 

N (T ) (±E ) N (h ) 1 (±E ) °^ ^ [multiplied by the factor i/(ui + iS)] are 

y ± = b( iw) + n-p( k+q) ^ b( _J Z k +q) m the paramagnetic and the diamagnetic term, respectively. 



iE n + h,uj =p i5 k 4 

In the preceding equations 



1E„ - JkJ =F -E^ 



1 



■S(q,w) = -Im{xsF(q,w)} 



(17) 



(18) 



is the spectral function of the spin fluctuations. The 
quasiparti cle excitation en ergies -E k are defined as usual: 
Ek = (e k — n) 2 + Aj, Nb and np are the Bosc- 
Einstein and Fcrmi-Dirac distribution functions, respec- 
tively. For a single Cu02 plane (1/iV) XZq m ^1 S - 13,16 
can be replaced with / 2DBZ a 2 d 2 q/(27r) 2 . 



B. Optical response 

The Kubo formula for the conductivity <Tjj(q, a;) (i, j 
e {a;, y, z}) reads 39 : 



(q, cj) = 



i<5 



Rij(ci,u)) H 4 



(19) 



where 



%(q,u;) = -^ / d*0(f)e tot <[fc(q,i)^(-q,0)]) 

— OO 

(20) 

is the retarded correlation function of the paramagnetic 
current density operator and no is the average electron 



The contribution to tr^ (q, uj) of charge carriers in a single 
band is given by a formula 40 which has the same structure 
as Eq. 19. The Fourier component of the corresponding 
current density operator is given by 



(21) 



where v(k) = (l/h)de/d\t, and the diamagnetic term is 
replaced with 



oj + 15 h 



2 V^ij ) ' Kij 



4e 



V ^ dhdk 



c+ a c ka . (22) 



For easy reference, we summarize here some basic 
properties of the optical conductivity <7y (to) = <xy(q = 
0,cj). It consists of a singular part a?j(u)), related to the 
condensate, and a regular part ^-(w). They are given 
by 



iS 



and 



a r (u) = - [R(u) ~ Re{i?(0)}] 



(23) 



(24) 



The tensor indices are omitted for simplicity, w p i iSC is the 
plasma frequency of the supcrfluid, 



w pl,sc — 



(K) - Re^O)} 

n 



(25) 
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and R(uj) = i?(q = 0, w). The imaginary part ^(w) of 
o{uj) can be expressed in terms of the real part, Ui(w), 
using the following Kramcrs-Kronig relation: 



i f 2 2 r°° 

<T2(?») = -~ 72 W + — P / 

ft 2 ' W7T J 0+ 



,u/ 2 CTi(a/) — W 2 CTi(w) 

(26) 

The symbol + in the lower limit means that the singular 
component of cr\, <J'{(oj) = keqW^ sc <5(o;) , is excluded from 
the integration. The function <Ji(uj) (including its singu- 
lar component) satisfies the following useful sum rule: 



Ir. 



du<7i(u) = -e uj 



pi 



Y ¥ {K), (27) 



which allows us to express the plasma frequency of the 
supcrfluid as 



^pl,sc 



2 



(28) 



7T£0 Jo+ 

The correlation function Rij{w), which is the key quan 



J 



tity of the theory, can be obtained by the analytical con- 
tinuation of the corresponding Matsubara Green's func- 
tion Pr/(q, ifl m )- Using the Wick's theorem and neglect- 
ing vertex corrections we obtain 41 : 



E v * ( k + V i ( k + L ( k ' * ifim), 



where i(k, q, i£l r 



(29) 



-— £ Tr[G(k + q,iE n + inn m )G(k,iE n )}. (30) 



n— — cxd 



The regular component of <xy i (w) can be expressed 
in terms of the matrix spectral function of Eq. 15 as 
follows 42 : 



Tijl 



|-i (k) Wj - (k) J ^Tr [A(k, S)A(k, E + hu)] x [n F (S) - n F (E + hu)] . 



(31) 



r 



superconductor containing N p equivalent 



For a 
weakly coupled Cu02 

(WEk in E q s - 22, 

(N p /d)f 2DBZ d 2 k/(2ir)i 
eter along the c-axis. 



and 

planes within a unit cell, 
29, 31 can be replaced with 
where d is the lattice param- 



III. INPUT PARAMETERS AND 
COMPUTATIONAL DETAILS 

A. Dispersion relation 

In this work we use a tight-binding expansion of the 
in-plane dispersion relation including the second nearest 
neighbor hopping terms: 

e(k) = — 2t[cos(k x a) + cos(k y a)] — 4t' cos(k x a) cos(k y a) . 

(32) 

The values of the parameters are the same as in Ref. 30 
and they are summarized in Table 1. The correspond- 
ing Fermi surface is shown in part (a) of Fig. 1. For 
comparison, the Fermi surface, corresponding to the 
six-parameter fit of the dispersion relation of optimally 
doped Bi2Sr 2 CaCu 2 8 (Bi-2212) used in Ref. 20, is 
shown in part (b). The figure is further used to introduce 
a notation of some important k-points. Within the SF- 
based scenarios, the quasiparticles around the CS (HS) 
are weakly (strongly) renormalizcd by the spin-fermion 



coupling. The CS coincides with the node of the super- 
conducting gap whereas the HS is located in the region 
of high values of | A | . 

Our approach will be applied here to the bilayer com- 
pounds, i.e., materials, that contain two Cu02 planes 
within a unit cell. The conduction band of these com- 
pounds can be expected to be splitted into the bonding 
branch and the antibonding branch 43 . For a long time, it 
has not been clear, whether such splitting indeed occurs 
or is supressed by strong correlations in the Cu02 planes. 
Two well defined bands were recently observed in photoc- 
mission experiments on strongly overdoped Bi-2212 and 
some signatures of the splitting have been reported also 
for optimally doped and underdoped Bi-22 1 2 44 - 45,46 (for 
a review see Ref. 21). It has been proposed that some of 
the effects associated earlier with the magnetic mode, in 
particular the peak-dip-hump structure of the photoemis- 
sion spectra, can be understood solely in terms of the two 
bands. Subsequently it has been shown, however, that a 
"true", i.e., linked to sclfencrgy effects, pcak-dip-hump 
structure develops in the spectra below T c 47 ' 48 . In the 
present study, the bonding-antibonding splitting is not 
taken into account, i.e., two identical conduction bands 
per unit cell are assumed, corresponding to two equiva- 
lent and independent planes. It is possible to go beyond 
this approximation 8,47,49 but we believe that it would not 
seriously modify our conclusions regarding the in-plane 
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a(A) d(A) 


t (eV) 


*' (eV) 


H (eV) A (eV) 


9 (eV) 


Tkjo (eV) 


r (eV) 


€(A) 


3.828 11.650 


0.250 


-0.100 


-0.350 0.030 


0.350 


0.040 


0.010/0 


9.0 (=2.35a) 



TABLE 1: Values of the parameters used in the computations. The values of a and d correspond to YBa2Cu307_<s. We 
consider two CuCh planes within a unit cell, as in YBa2Cu307_^, i.e., N p = 2. 
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FIG. 1: Fermi surface (thick line) for the dispersion given by 
Eq. 32 (a) and for the dispersion of Ref. 20 (b). Supplemen- 
tary contours (thin lines) correspond to excitation energies of 
±0.033 eV and ±0.050 eV. A notation of some special points is 
introduced. The point CS ("cold spot") is the crossing point 
of the Fermi surface and the BZ diagonal (dashed-dotted line 
starting in the T point). The point HS ("hot spot") is the 
crossing point of the Fermi surface and the antiferromagnetic 
BZ boundary (dashed-dotted line starting in the X point). 
The point B is located between the CS and the HS and its 
(dimensionless) coordinates are (1.96, 0.60) and (1.96,0.78) in 
(a) and (b), respectively. 



infrared spectra. 



B. Superconducting gap 

If not specified in the text, we use a common ansatz 
for the superconducting gap of d a 2_ y 2-symmetry: 

A(k) = -^-[cos(k x a) — cos(fcj,a)] (33) 

with Ao = 30 meV, a value close to the experimental 
result for optimally doped Bi-2212 21 . The coupling to 
the spin-fluctuations slightly enhances the value of the 
gap in the region around the HS. For the value of g 
of 350 meV used in our computations, the renormalized 
value of | A | exactly at the hot spot is about 35 meV. This 
fact, however, has no significant impact on the quantities 
discussed. 



C. Spin susceptibility and the coupling constant 

In most of the computations we use the spin suscepti- 
bility of the resonance-mode form of Eq. 4 with the values 



of the input parameters given in Table 1. The FWHM 
of the resonance peak T is small as compared to both loq 
and A . It can be therefore neglected, in the first approx- 
imation, an approach applied by Eschrig and Norman 20 . 
We have performed the calculations for both T = 10 meV 
and T = 0, i.e., 



Fit 



i±(q-Qre 2 



6(u> - uj ) , 



(34) 



in order to identify the role of the finite linewidth. In 
the following we shall refer to the calculations with F = 
lOmeV and r = as to those with the broad and the 
sharp magnetic resonance, respectively. The value of the 
"oscillator-strength" parameter F is determined by using 
the same normalization condition as in Rcf. 30: 



hi = h 



i f$ r dLj (2 NB{hui) + 1} jB(q ' w) = 2 



(35) 

where B(q,ui) is the spectral function of the spin fluc- 
tuations defined by Eq. 18. Equation 35 represents an 
approximate "total momentum sum rule" (see Ref. 18 
for a discussion). 

Note that when using xsf of the resonance-mode form 
together with Eq. 35 we make two approximations, (a) 
Within the model of two identical conduction bands per 
bilayer unit cell the effective spin susceptibility is equal 
to the average of the so called odd (xo) and even (xe) 
component (for a definition see Ref. 18). Here we take 
into account only the resonant part of Xo, i.e., we do not 
consider % e and the continuum part of Xo- This may be a 
reasonable starting point for an analysis of the low-energy 
phenomena. First, the magnitude of Xo is typically by a 
factor of 2-3 larger than that of x e 18 ' 50 , and second, the 
continuum part cannot be expected to cause any sharp 
structure in the quasiparticle selfenergy (see Ref. 20 for 
a detailed discussion), (b) In Eq. 35 the resonance is 
assumed to collect one half of the total magnetic spectral 
weight, Im = 1/4. Experimentally, the value of Im in 
Eq. 35 corresponding to the mode is only about 0.01 x 
(1/4) for optimally doped YBa 2 Cu 3 7 _ 5 (Y-123) 18 and 
only ca 0.04 x (1/4) for optimally doped Bi-2212 17 . Since 
F ~ Im enters Eq. 2 in the product g 2 F, however, we 
continue to use Eq. 35 for simplicity keeping in mind 
that the actual value of the coupling constant required is 
an order of magnitude higher than the one of Table 1. 

Let us make here a short comment on a recent discus- 
sion regarding the small experimental values of Im ■ Kee, 
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Kivelson, and Aeppli 51 claimed that they are too small 
for the mode to play any significant role. Abanov et al. 
in their response to the criticism argued that even though 
Im is small, the mode can lead to large selfenergy effects 
because it is rather narrow in q space. There is another 
argument justifying the "mode models" . As mentioned 
above F ~ Im appears in Eq. 2 only in the product g 2 F. 
The susceptibility xrm normalized to 1/4, as in our cal- 
culations, together with g = 0.35 cV, will be shown to 
lead to results that are to some extent in agreement with 
experimental data. The same results could also be ob- 
tained using "the true xrm", i-e., Xrm normalized to 
the experimental value of 0.002 (0.009), with g sa 4eV 
(g f=a 2 eV) . The important point is that the latter values 
of g are not unrealistic given Eq. 1 and the high values of 
U expected to be relevant for the HTCS. It is thus well 
possible that, inspitc of the small values of Im, the selfen- 
ergy effects due to the mode are large. The fact, however, 
that its spectral weight is small as compared to that of 
the SF continuum implies that only some low energy fea- 
tures can be attributed to the mode. The intermediate 
energy-scale spectra must be determined rather by the 
continuum. 



D. Computational details 

The results presented in this work have been obtained 
using two different codes. In the first one, the finite 
linewidth of the resonance is taken into account, the sec- 
ond one is adapted to the case of the sharp resonance. 
The calculations proceed essentially in three steps. 

(i) In the first step, the self-energies and spectral func- 
tions of selected k-points are computed by using Eqs. 16 
and 15. First we concentrate on the frequency integra- 
tion of Eq. 16. Within the sharp resonance approach, 
it is trivial due to the S function in Eq. 34. The struc- 
tures of the resulting function, however, are sharp and 
the sampling of the subsequent BZ integration must be 
very dense. Within the broad resonance calculation, the 
integral of Eq. 16 is divided into two parts: the part con- 
taining the Bosc factors in the numerators of Eq. 17 and 
the one containing the remaining terms. The latter part 
can be expressed analytically but the former one has to 
be calculated numerically. For this purpose, we have used 
the grid containing 500 ^-points ranging from lu = to 
u> = 10kBT/h. Due to the separable form of xrm, the 
same frequency integrals occur for any k g l.BZ. This 
has allowed us to compute their values once and for all 

I 



at the beginning. 

Next, we focus on the BZ integration. The dispersion 
relation is nontrivial and the structure in the suscepti- 
bility relatively sharp. In addition, a discussion of the 
singularities requires the spectral functions to be contin- 
ued very close to the real axis: a typical value of S in 
Eq. 15 is only 0.05 meV. For all these reasons a very 
fine grid (typically 300 x 300 for the broad resonance and 
1000 x 1000 for the sharp one) is needed for the BZ inte- 
gration. Without further simplifications the calculations 
of the infrared conductivity would be fairly time con- 
suming: they would require elaborate integrations for all 
k-points of the corresponding grid. This problem can be 
partially overcome by using an approximation described 
below. The BZ integral of Eq. 16 can be schematically 
written as 

/ dq£(q)I?(k + q,^„), (36) 
J 2D BZ 

where £ is the q dependent part of B and D is a function 
of k + q = k' and iE n . Using the k-space periodicity of 
both £ and D we obtain 

/ = / dk'C(k'-k) J D(k',^„). (37) 

J 2D BZ J 2D BZ 

The idea underlying our approximation is that £ is a 
smooth function as compared to D. This allows us to 
use the approximation 

/ «^£(]%-k)Ai(i£y, (38) 

J 2D BZ tj 

where i = 1,2, ...Ni, j = 1,2, ...Ni, N\ specifies the 
number of k-points in the auxiliary grid (typically Ni = 
30), ky = (kxukyi), k xi = (2i - 1)tt /(aN{), k yj = (2j - 
l)it/(aNi), and 

Dij(iE n ) = / <5kD(k,iE n ). 

J2(i-l)ir/(aNi) J2(j-l)n/(aN 1 ) 

(39) 

The quantities Dij(iE n ) have been computed once and 
for all at the beginning, using grids of typically 10 x 10 
(30 x 30) k-points for the broad (sharp) resonance. 

(ii) In the second step, the contributions of selected k- 
points to the real part of the (isotropic) in-plane infrared 
conductivity are calculated. We define the contribution 
of a k-point as (cf. Eq. 31) 



Aai(k, W ) = |-^[t£(k)+i;2(k)] / ^Tr[A(k,E)A(k,E + hcj)] x [n F (E) - n F (E + M] • (40) 



The computation of the convolution of the spectral func- 



tions in Eq. 40 is again very demanding because the sharp 
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quasiparticlc peaks must be handled correctly. Typically 
a grid of at least 5000 energy-points for a spectrum rang- 
ing up to 300 mcV is required. In the calculations of 
Subsecs. IV E and IV F, where the sharp spectral struc- 
tures are not at the centre of interest, the spectral func- 
tions have been continued only to a distance of typically 
1 mcV from the real axis and energy grids of typically 
6000 points for spectra ranging up to 3 eV have been 
used. 

(hi) The calculation of cf\{lS) is completed by the BZ 
integration of Acti : 



<ri(w) = 4a 2 



dkAai(k,a;) 



(41) 



2DIBZ 



where IBZ stands for the irreducible part of the BZ. The 
integration is relatively simple because the functions Acri 
are already fairly smooth. The grid of 200 x 200 k-points 
in the full BZ, i.e., ca 5000 k-points in its irreducible 
part, proved to be sufficient. 

If not specified in the text, the computations are per- 
formed for T = 20 K. 



IV. 



RESULTS AND DISCUSSION 



The infrared spectra will be analyzed in terms of the 
quasiparticle spectral functions. For this reason, we first 
present and discuss, as a necessary prerequisite, our cal- 
culated spectral functions of some representative k-points 
(Subsec. IV A) and their contributions to the real part of 
the conductivity (Subsec. IV B). The main new results of 
the paper will then be presented in Subsects. IVC-IVF. 



A. Quasiparticle spectral functions 

Figure 2 shows the spectral functions (more precisely, 
the first diagonal elements of the matrix (15), i.e., the 
"true" spectral functions) of three points located at the 
Fermi surface (the points CS, B, and HS) and of the X- 
point. The spectra exhibit quasiparticle pcak(s) (one for 
the CS and two for the other points, where the magni- 
tude of the superconducting gap is finite) , that would ap- 
pear already for the (noninteracting) BCS ground-state, 
and incoherent broad satellites at higher energies, that 
are due to the interaction. The structures are more pro- 
nounced in the spectra for T = 0. The latter are thus 
more suitable for a discussion of the singularities. 

The energies of the quasiparticle peaks labelled Q and 
Q' (see Fig. 2) are — 8meV (CS, the nonzero value is due 
to a slight renormalization of the Fermi surface caused 
by the interaction), ±22meV (B), ±35meV (HS), and 
±37meV (X). The magnitude of the values for the HS is 
larger than Ao which reflects an enhancement of super- 
conductivity due to the coupling to the mode. 

Next we discuss the incoherent parts of the spectral 
functions. For all of the k-points, three onset features can 
be resolved at negative energies, denoted by R, S, and T 
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CS 



r = o 

r = 1 meV 
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FIG. 2: Spectral functions of the k-points defined in Fig. 1. 
The notation of spectral structures is introduced in panel (b). 



(see Fig. 2), and two at positive energies (R', S'). Interest- 
ingly, their positions are the same for all the k-points: the 
feature R(R') is located at -(+)hu [T-40meV], S(S') at 
-(+)(fiwo + A A ) ~ -(+)(Rwo + A ) [« T-70meV] and T 



at -\hu + E x \, where E x = y/e x + A 2 X [w -lOOmeVj. 
The physical meaning of the structures S,S', and T has 
been discussed in detail by Eschrig and Norman 20 ; their 
analysis also yields the formulas for the energies of the 
features. The structures R and R', which are well re- 
solved only for T = 0, correspond to final states involving 
a nodal quasiparticle and the magnetic excitation. The 
energy of the broad maximum at negative energies, below 
the feature T, is k-dependent. Its values are — 100 me V 
(CS), -115 meV (B), -150 mcV (HS), and -180 meV 
(X) . The decrease of the energy when going from the CS 
to the X-point is accompanied by a systematic increase 
of the spectral weight. 

In order to illustrate the role of the input parameters, 
we compare in Fig. 3 the results presented above with 
those obtained using the parametrizations of the disper- 
sion relation and of the magnetic susceptibility and the 
values of A (35 mcV), g = 0.65, and T (40 K) of Ref. 20. 
Note first that the reference spectra are in good agree- 
ment with those of Ref. 20. Note further, that for small 
absolute values of energy, the two sets of spectra are very 
similar. On the other hand, the positions and shapes of 
the lowest incoherent peak are fairly different. This is 
due to the fact that the value of the effective coupling 
constant used by Eschrig and Norman is considerably 
smaller than the ours 53 . The difference in the energy of 
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the onset of the latter peak (ca — 100 me V in this work, 
ca —90 meV in Ref. 20) is related to the difference in the 
location of the X-point with respect to the Fermi surface 
(see Fig. 1). 
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FIG. 4: Contributions of the selected k-points to the real part 
of the optical conductivity. The onset features are marked by 
the arrows. 



FIG. 3: Spectral functions from Fig. 2 (solid lines) compared 
with those obtained using the values of the input parameters 
of Ref. 20 (dashed lines). 



B. Contributions of the selected k-points to a\(uS) 

The contributions of the k-points labelled as CS, B, 
and HS to a\{ui) are shown in Fig. 4. The contribution of 
the X-point is identically zero because v(X) = 0. In the 
following, the structures of the functions will be discussed 
in terms of those of the corresponding spectral functions. 
Such a discussion is straightforward only for the normal 
state, where the expression on the r. h. s. of Eq. 31 re- 
duces to a convolution of the (true) spectral functions. 
For the superconducting state, the off-diagonal compo- 
nents of the spectral-function matrix (15) play a cru- 
cial role. For example, the transitions Q — * Q' do 
not contribute because of a precise cancellation of the 
term on the r. h. s. of Eq. 31 containing the diagonal 
components and that containing the off-diagonal com- 
ponents. Other transitions between occupied and unoc- 
cupied states, however, do show up in the conductivity 
spectra, only their spectral weights are influenced by the 
presence of the off-diagonal components. 

Three singularities (discontinuities of the first deriva- 
tive), can be observed in the spectra. In Fig. 4 they 



are marked by the arrows. The onset feature Oi corre- 
sponds to the transitions Q — > R' and R — * Q' (for the 
CS, Q = Q'). It appears at the energy of + huJo, 
where E^ is the energy of the quasiparticlc peak. The 
values of the corresponding wavenumbers are 380 cm , 
500 cm -1 , and 600 enr 1 for the points CS, B, and HS, re- 
spectively. The second onset feature O2 is related to the 
transitions Q — > S' and S — > Q'. It is shifted with respect 
to 0\ approximately by An, the values of the correspond- 
ing wavenumbers are 620 cm -1 , 740 cm -1 , and 840 cm -1 , 
respectively. The third discontinuity of the first deriva- 
tive labelled O3 can be traced back to the transition 
T — > Q'. It is shifted with respect to Oi by Ex, the 
values of the corresponding wavenumbers are 850 cm -1 , 
970 cm -1 , and 1070 cm -1 , respectively. Small structures 
at very low energies in parts (b) and (c) of Fig. 4 are 
artefacts due to an incomplete continuation to the real 
axis. The contribution of the cold spot requires a more 
careful discussion. Besides the features Oi, O2, and O3, 
it contains a Drudc-likc component at very low energies. 
Its presence is related to the fact that at this point the 
superconducting gap vanishes. Some structures around 
60 cm -1 and 320 cm" 1 are due to transitions between the 
maximum of the quasiparticle peak Q located at —8 meV 
and its tail ranging to E = 0, and between the tail and 
the structures R,R', respectively. The spectral weight as- 
sociated with the transition T — > Q' is small because only 
the tail of the quasiparticlc peak participates. It can be 
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seen in Fig. 4 that the spectra for finite T are similar to 
those of r = 0, the onset features, however, are by far 
less pronounced. 
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FIG. 5: Contributions of the selected k-points to the real part 
of the optical conductivity from Fig. 4 (solid lines) compared 
with those obtained using the values of the input parameters 
of Ref. 20 (dashed lines). 

For reference, we compare in Fig. 5 the results pre- 
sented above with those obtained using the values of the 
input parameters of Ref. 20 (the values of a, d and 6 have 
not been changed) . It can be seen that the frequencies of 
the singularities are fairly close to each other, except for 
O3. The discrepancy can be traced back to the difference 
in the location of the Fermi surface. The differences in 
the absolute values reflect those of the dispersion relation 
and of the coupling constant. 

The contribution of the hot spot shown in the bottom 
panels of Figs. 4 and 5 can be compared with <ji(lo) re- 
sulting from the model of ACS, shown in Fig. 2 of Ref. 33. 
The latter exhibits a sharp onset feature at 2Ao + hujQ, a 
maximum around 5 Ao , and a gradual decrease at higher 
frequencies. Our spectra display all these features. The 
onset is labelled as 02- The maximum is located at about 
1500 cm" 1 for our values of the input parameters and at 
ca 1100 cm -1 for those of Ref. 20, its frequency strongly 
depends on the value of g. The main differences are: (a) 
In the spectra of Ref. 33, 01 = for frequencies below 
the onset. In our spectra, Aa± acquires nonzero values 
already above the first onset frequency of ~ Ao + Tlloq. 
The discrepancy is due to the fact that ACS neglected 
excitations involving nodal quasiparticles. (b) In our 
spectra, the onset O2 is followed by another one at ca 



Ao + Ex + huo (O3) . Its absence in the spectra of Ref. 33 
is due to the fact that details of the dispersion relation 
have been neglected. 

In summary, the contribution of a k-point to o\ (uj) 
exhibits three onset features, whose energies are E^+fiuo, 
ca + Ao + frwo, and E^ + Hluq + Ex- 



C. Real part of the optical conductivity and the 
origin of the onset feature around wo 

So far we have discussed singularities that appear in 
the contributions to of the individual k-points. It 

is important to find out how do these singularities man- 
ifest themselves in the spectra of the total conductivity. 
Here we shall show that the feature Oi of Fig. 4 gives 
rise to a clear onset of o\ starting around too- The singu- 
larities labelled as O2 and O3 cannot be clearly resolved 
in the spectra themselves, the former one, however, can 
be uncovered by means of numerical differentiation. This 
will be shown in the next subsection. 

Figure 6 shows the spectra of <J\(u) (upper panel) and 
the contributions of the three regions of the irreducible 
BZ defined in Fig. 7 (bottom panel). The imaginary part 
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FIG. 6: Real and imaginary parts of the optical conductivity 
(upper panel) and the contributions to ai(iv) of the three 
regions of the irreducible BZ defined in Fig. 7 (bottom panel). 

(T 2 (w) has been obtained by using Eq. 26. The values of 
Wpi and Wpi, sc obtained by using Eq. 27 and Eq. 28 are 
17800 cm -1 and 12100 cm -1 , respectively Interestingly, 
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FIG. 7: Three regions of the irreducible BZ used in the dis- 
cussion of the total conductivity. 
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the magnitudes of the three contributions to ai(u>) in the 
infrared are comparable. This is a result of a competi- 
tion between the velocity factor in Eqs. 31, 40, which is 
higher for the quasiparticles around the BZ diagonal (re- 
gion III), and the density of states at the Fermi surface, 
which is higher in region I. It appears that the charac- 
teristic singularities discussed in the previous subsection 
can still be resolved in the spectra of the bottom panel. 
For example, the contribution of region I still contains 
the onset features Oi, O2, and O3 characteristic of the 
hot spot. In the spectra of the total conductivity, how- 
ever, they are approximately averaged out and cannot be 
clearly resolved, except for the first onset feature: a\ is 
rather small above the narrow Drude peak and below u>q 
and it starts to increase at lo (or, for finite T, around 
ujq). The bottom panel of Fig. 6 shows that the onset 
originates in region III. It can be attributed to the lowest 
allowed transitions (Q — ► R' and R — ► Q') of the region 
around the nodes. The shape of the onset can thus be 
expected to depend on the behavior of Ak in the proxim- 
ity of the nodes. In order to explore this dependence, we 
have performed the computations for the following three 
forms of the superconducting gap: 

A k = ^[cos^a) - cos(fc. y a)] 7 (42) 

with 7 = 0.5, 1, 3 (for 7 = 0.5 the signs have to be prop- 
erly adjusted). The results are shown in Fig. 8. They 
can be interpreted as follows: the larger (smaller) the 
value of the exponent 7, the wider (narrower) the region 
around the BZ diagonal, where the magnitude of the su- 
perconducting gap is small, and the larger (smaller) the 
spectral weight right above loq originating in this region. 
The increase (decrease) of 7 thus effectively results in a 
red (blue) shift of the onset. 

For reference, we again compare in Fig. 9 the results 
presented above with those obtained using the values of 
the input parameters of Rcf. 20; the corresponding values 



of w p i and ui p ^ : 
tively. 



are 17900 cm and 12700 cm , respec- 



D. Parameters of the extended Drude model, 
effective spectral function and the origin of its 

main maximum 

A common way to discuss the in-plane infrared spec- 
tra is in terms of the scattering rate [1/t](w) and the 
mass enhancement factor [m* /m](u>). These quantities 
arc defined by the so called extended Drude model 24,32 : 



[1/t](uj) — iuj[m* /m](uj) 



(43) 



Their physical significance relies on a close relation be- 
tween the optical selfenergy, 



£ op *(w) = 



1 



1 



-M 



--M 

T 



(44) 



and the (true) quasiparticlc selfenergy (for a recent dis- 
cussion see Ref. 54). The calculated spectra of 1/r and 
in* jm are shown in Fig. 10. They have been extracted 
from those of o\ and 02 using the following formulas re- 
sulting from Eq. 43: 



T 



-M 



2 TlM 



(72 (u) 



(45) 



(46) 



The trends in the spectra of 1/r and m* jm can be eas- 
ily related to those of <7i and 02 . The obvious similarity 
between the spectra of 1/r and those of 01 (w), for exam- 
ple, can be interpreted as follows. At low frequencies (a) 
(T2 3> "l and (b) 02 is dominated by the singular term, 
£ o w pi,sc/ w ( see Fig- 6). As a consequence, [1/t](w) is ap- 
proximately proportional to lo 2 <ji(u)) (see Eq. 45). This 
explains the similarity. The latter allows us to consider 
derivatives of [l/r](w) instead of those of (J\(ui) when 
searching for the singularities. 
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FIG. 10: Spectra of the scattering rate and the mass enhance- 
ment factor extracted from those of o\ and <j% shown in Fig. 6. 



Recently, it has become popular 



26.31.55,56,57,58 



to use 



the function W(u>) of Eq. 6 which is related to [1/r] as 



W{ui) 



27rdw 2 



u) — (u)) 

T 



(47) 



This has been motivated by the theoretical works of 
CSB 31 and ACS 33 discussed in the introduction. We re- 
call that according to CSB the most pronounced maxi- 



mum of W(u>) reveals the spectral function of the bosons, 
that are coupled to the charged quasiparticles (shifted by 
Ao, see Eq. 6). Within the resonant mode scenario, the 
maximum should thus be located at f2csB = + Two- 
According to the strong-coupling SF-based theory of ACS 
W(u>) should exhibit three singularities at 2Aq + huio, 
4A , and 2A + 2hcj (see Fig. 3 of Ref. 33). The first 
singularity, which is the most pronounced one, is due to 
the first onset of ai(u>) (see Fig. 2 of Ref. 33 and the 
discussion in Subsect. IV B). It consists of a sharp max- 
imum followed by a minimum. 

Here we shall follow the strategy of the preceding sub- 
sections and attempt to interprete the structures of W(uj) 
in terms of the quasiparticle spectral functions and the 
corresponding contributions to a\{ui). Our most impor- 
tant finding is that the main maximum of W(ui) can be 
associated with the onset feature O2 of parts (a) and (b) 
of Fig. 4, i.e., with the appearence (above the character- 
istic frequency of the maximum) of final states consisting 
of a quasiparticle from the nodal region, a quasiparticle 
from the antinodal region, and the resonance mode. 

Figure 11 shows our calculated spectra of W{ui) (up- 
per panel) and contributions to W of the three regions 
defined in Fig. 7 (bottom panel). The contribution of a 
region is defined by 



AW(lu) 



£o^pi d 2 
2tt du 2 



A<7l(w) 



(48) 



The spectra of the quantities in the square brack- 
ets of Eqs. 47 and 48 have been smoothed before 
taking the second derivatives, a typical scale of the 
smoothing being 3meV. First we discuss the spectra 
of the bottom panel. Each contribution exhibits three 
maxima (Mi, M 2 , M 3 ) separated by two minima (M{, 
M' 2 ). The corresponding frequencies are 800 cm -1 , 
1090 cm" 1 , 1800 cm" 1 , 910 cm" 1 , 1500 cm" 1 (I), 
680 cm" 1 , 950 cm" 1 , 1620 cm" 1 , 830 cm" 1 , 1130 cm" 1 
(II), 600 cm" 1 , 880 cm" 1 , 1030 cm" 1 , 770 cm" 1 , 
940cm" 1 (III). Spectral features above the frequency 
of the third maximum are weak. It can be relatively 
easily seen that the maxima M\ and M2 correspond 
to the onset features O2 and O3, respectively (see 
Figs. 4 and 6) and that the mimimum M' 2 is related to 
the maximum of A<7i(w) in mid-infrared (see Fig. 6). 
The structures M[ and M 3 can be interpreted analo- 
gously. Next we concentrate on the upper panel. The 
spectra of W(u>) exhibit a pronounced maximum at 
~ 620 cm" 1 , a minimum at ~ 870 cm" 1 , weak structures 
in the range from 1000 cm" 1 to 1400 cm" 1 , and a 
minimum around 1500 cm" 1 . All these features can 
be interpreted in terms of constructive or destructive 
interference of the three contributions. In particular, the 
main maximum originates from the structures Mi (III) 
and Mi (II), the latter being partially compensated 
by M{ (III). Note the important role played by the 
negative contribution of the zone III. The frequency 
of the maximum CIm (here 620 cm" 1 ) appears to be 
slightly higher than w(Mi,77/) w u(0 2 ,UI) (here 
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FIG. 11: Spectra of the function W(uj) defined by Eq. 6 ob- 
tained from the conductivities of Fig. 6 (upper panel). Bot- 
tom panel: the contributions of the three regions defined in 
Fig. 7. 



600cm : ). The latter frequency is further slightly 
higher than cj(0 2 ,CS) w (A /h) + iv (565cm- 1 ). 
Our general conclusion is that Qm is somewhat larger 
than SIcsb = (Aq/K) + but much smaller than 
2(Aq/K) + u>o- The magnitude of the difference, 
fifcf — (Ao/h) — uiq, depends on details, e.g., on the 
shape of the superconducting gap. In our case we have 
Qm — (Ao/h) — ljq = 55 cm -1 . The conclusion should 
be valid for any sharp bosonic mode which couples the 
nodal regions to the antinodal ones. It can be easily seen 
that the shape of the maximum does not correspond to 
the neutron peak. The FWHM is about 170 cm- 1 and 
about 210cm- 1 for T = and T = lOmcV (ca 80cm- 1 ), 
respectively. 

Figure 12 shows the results obtained using the values 
of the input parameters of Ref. 20. The main differences 
between the spectra of the bottom panels of Fig. 11 and 
Fig. 12 are: (a) The maxima Mi and are closer to 
each other in Fig. 12 than in Fig. 11. (b) The maxi- 
mum M3 is located at lower energies in Fig. 12 than in 
Fig. 11. These differences can be traced back to those 
in the position of the Fermi level and in the value of the 
coupling constant. The frequency of the main maximum 
in the upper panel of Fig. 12 is approximately the same 
as that of Fig. 11, in agreement with our general conclu- 
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FIG. 12: (a) The T = OmeV spectrum of W(uj) from Fig. 11 
(solid line) together with the one corresponding to the val- 
ues of the input parameters from Ref. 20 (dashed line), (b) 
The same as in Fig. 11 (b) but for the values of the input 
parameters from Ref. 20. 



sion presented above. The shape of the maximum and 
the spectral structures at higher energies, however, differ 
considerably. 



E. Role of the spin-fluctuation continuum 

So far we have considered the spin susceptibility in the 
resonant-mode form of Eq. 4. Such an approach may 
suffice to provide an understanding of the spectral fea- 
tures discussed in the preceding subsections: the onset 
of <J\(ui) around ojq and the main maximum of W(ui). It 
is well known, however, that the spectral weight of the 
resonant mode is much smaller than that of the spin- 
fluctuation continuum. The overall shape of the opti- 
cal response, in particular in the mid-infrared region, 
can thus be expected to be determined by the contin- 
uum rather than by the mode. The inadequacy of the 
rcsonant-modc approach to reproduce the mid-infrared 
data can indeed be seen by comparing the computed scat- 
tering rate spectra with those deduced from experimen- 
tal data 27 ' 29 and called "experimental (scattering rate) 
spectra" in the following 59 . In the spectra of Fig. 10, 1/r 
saturates above 1200 cm -1 at a value of about 1200 cm -1 
and it decreases at higher frequencies (to be shown be- 
low) . In the experimental spectra, on the other hand, 1/r 
saturates only above ~ 5000 cm -1 , at a value of several 
thousands cm . 

Here we report on our attempts to understand the ex- 
perimental data in terms of a model where the spin sus- 
ceptibility contains, in addition to the resonant mode, a 
very broad continuum term: 



X(q, w) = &MXRM (q, w) + b c \c (q, . 



(49) 
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where XBM is given by Eq. 4 with T 
junction with Eq. 35 and 



xc(q,w) 



1 



i + (q-Q) 2 e 2 



lOmcV in con- 



(50) 



c 



iT c uj 



is the continuum component with the same structure as 
Xrm- The values of the parameters of xc used in our cal- 
culations arc: £c = 2 A (i.e., about one half of the value 
of the lattice parameter), u>c = 400 me V (this choice will 
be motivated below), = 1000 meV, and the value of 
Fq is determined by Eq. 35. The spectra of both com- 
ponents are shown in Fig. 13 (a). A similar continuum 
component has been considered by Schachinger, Tu, and 
Carbotte 37 in their discussion of the scattering rate in 
Bi-2212. The parameters 6m and 6c in Eq. 49 express 
the spectral weights of the two components. Figure 14 



■ 

03 



(a) 



, 1 1 , 1 1 , r- 

resonance mode 
continuum (10 x) 



1 1 1 1 1— 



(b) 



SC (T=20 K) 
N(T=100K) 



200 



400 600 
E (meV) 



800 



1000 



FIG. 13: (a) The (q-integrated) spectra of the imaginary part 
\" of the spin susceptibility corresponding to the resonance 
mode (solid line) and to the continuum term introduced in 
Eq. 50 (dashed line), respectively, (b) Those used as a starting 
point in the computations of Subsec. IV F. The solid and the 
dashed line correspond to the superconducting state (T — 
20 K) and to the normal state (T = 100 K), respectively. 

shows the 6m- and 6c-dependence of ai(ui) and [l/r](w). 
The low frequency parts of the spectra shown as the solid 
lines in panels (ai) and (&2), have already been presented 
above in Figs. 6 and 10. Note that in the mid-infrared 
region (a) a\ decreases with increasing frequency much 
faster than in the experimental spectra and (b) 1/r also 
decreases, in contrast to the experimental spectra. The 
right panel of Fig. 14 clearly demonstrates the distinct 
roles of the two components: (i) The mode yields a rela- 
tively sharp onset of 1 /r centered at about f2jyr , followed 
by a gradual decrease, (ii) The continuum leads to an ap- 
proximately linear increase followed by a plateau above 
~ 4000 cm -1 . As shown in Fig. 15, the onset frequency 




FIG. 14: The conductivity and the scattering rate as functions 
of the parameters 6m and foe expressing the spectral weights 
of the mode and the continuum, respectively. 



of the plateau reflects the width of xc , which is here de- 
termined by loq of Eq. 50. Our choice of the value of 



c 

£■ 

CO 
CO 



co^OO meV 
co^OO meV 
co^OO meV 
co^OO meV 




(b) 




- / / 




i i ■' . 

■ / / / 

if// 








/ / // 




///. 





1500 E 



500 
E (meV) 



4000 
V (cm" 1 ) 





8000 



FIG. 15: The (q-integrated) imaginary part x" °f the spin 
susceptibility (a) and the scattering rate (b) as functions of 
the parameter u)c introduced in the text. 

ojq has been motivated by the experimental spectra of 
Refs. 27 and 29. (iii) It can be seen that 1/r is approx- 
imately a linear function of 6m and be- For example, 
[1/r] (6a/ = 1,6c = 4) [dotted line in part (02)] is ap- 
proximately equal to [1 /r] (6m = 1, be = 0) [solid line in 
part (a 2 )] + [1/t](6j\./ = 0, be = 4) [dashed line in part 
(ba)]. The value of 6a/ (6c) determines the value of 1/r 
above the sharp onset or, in other words, the height of 
the step in 1/r (the difference between the value of 1/r 
of the plateau and the height of the step) . (iv) The spec- 
tra for 6a/ 7^ exhibit an additional structure around 
2500 cm -1 . It can be shown, along the lines of Sub- 
sect. IV B, to be connected to a singularity in the density 
of the relevant final states located approximately at 



[i - e r + A + hu) . 



(51) 
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The latter is related to the van Hove singularity associ- 
ated with the bottom of the band at the r-point. The 
frequency us can thus be used to estimate fj, — er, i.e., 
the width of the occupied part of the band. It is very 
interesting that the experimental scattering rate spectra, 
both for Y-123 and for Bi-2212, display a fine structure 
around 3000 cm -1 , which is similar to that appearing in 
the computed spectra. The similarity leads us to spec- 
ulate that it has the same origin, i.e., that it is due to 
the singularity of the density of states discussed above. 
The value of fi — er of ~ 0.3 eV resulting from Eq. 51 is in 
rough agreement with the photomession values of 0.25 eV 
(Y-123, Ref. 60) and 0.40 eV (Bi-2212, Ref. 21). Accord- 
ing to the proposed interpretation, the structure should 
become less pronounced above T c because of the disap- 
pearance of the sharp mode. In the experimental spectra, 
however, it seems to be only weakly temperature depen- 
dent. At present we are not aware of any explanation of 
this discrepancy. 

It can be seen that the experimental spectra can be 
approximately reproduced only by taking a value of 6m 
of around 1.0 and a considerably higher value of 6c- For 
6m = 1.0 and 6c = 4.0, e.g., the calculated spectra are 
in reasonable agreement with those of optimally doped 
Y-123 from Ref. 29. In order to reproduce the spectra 
of optimally doped Bi-2212 from Ref. 27 a higher value 
of 6c (ca 10) has to be taken. Let us note that the 
values of the ratio 6c/6m under consideration are still 
smaller than those suggested by the neutron scattering 
experiments 17,18 . 
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F. Some aspects of the temperature dependence of 
the spectra 

In Subsubsect. IV Fl we present and discuss a com- 
puted temperature dependence of the spectra character- 
istic of the spin-fluctuation scenario, which may be help- 
ful for understanding the experimental data. In Subsub- 
sect. IV F 2 we then focus on the temperature dependence 
of the intraband spectral weight. 



1. Temperature dependence of the infrared spectra 
characteristic of the spin-fluctuation scenario 

Figure 16 shows the spectra of <7i, 1/r and m*/m com- 
puted for three temperatures: T = 20 K (superconduct- 
ing state), T = 100 K, and T = 200 K (normal state). 
For each temperature, the value of the chemical poten- 
tial was adjusted to yield the total number of electrons 
of 0.76 per unit cell: \i = — 355.5 mcV, — 347.5 mcV, and 
-345.7 meV for T = 20 K, 100 K, and 200 K, respectively. 
The values of the other input parameters for the nor- 
mal state are the same as those for the superconducting 
state given in Table 1, except for Ao = 0, T = 70meV 
(the mode is assumed to be overdamped), and £ = 6 A. 
Both for the superconducting and for the normal state 



FIG. 16: Calculated temperature dependencies of the conduc- 
tivity, the scattering rate, and the mass enhancement factor. 
The dashed-dotted line in part (b) corresponds to an unre- 
alistic case of the normal state, T — 100 K, and the spin 
susceptibility containing the sharp mode. 

we furher include the continuum contribution to x i n the- 
form introduced in Subsect. IV E, with four times larger 
spectral weight than that of the mode (6c /6m = 4). The 
20 K and 100 K spectra of \" are shown in Fig. 13 (b). 
The calculated values of w p i and cj p i )SC are 17500 cm -1 
and 10500 cm -1 , respectively. 

It can be seen that the model is capable of reproducing 
not only the main characteristics of the low temperature 
infrared spectra but also, at least qualitatively, the differ- 
ences between the normal and the superconducting state. 
Our discussion will be limited to the following issues: 

(i) The frequency dependence of Aa± = (Ti(100K) — 
(Ti(20K) in mid-infrared. In both the experimental (sec, 
e.g., Ref. 29) and the computed spectra (see Fig. 16), Aoi 
is a positive and (approximately) decreasing function of 
frequency. In the experimental spectra the decrease is 
very rapid below ~ 1000cm -1 , i.e., in the region of the 
spectral gap. Above 1000 cm -1 , the superconducting and 
the normal state spectra are close to each other and Aeri 
is very small. In the computed spectra, on the other 
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hand, the decrease of Aeri is relatively gradual. The dis- 
crepancy is very probably due to the fact that our ap- 
proach is not self-consistent. We believe that it would 
vanish or diminish if a sclfconsistent version of Eq. 13 
(i.e., Go replaced with G) was used. 

(ii) The onset of 1/r. The 100 K spectrum does 
not display the sharp onset around 700 cm -1 discussed 
above. Nevertheless, at low frequencies it deviates 
slightly from the straight line that can be used to ap- 
proximate [1/t]((x>) in the region between ca 1000 cm -1 
and 2500 cm" 1 , an effect similar to the one of some exper- 
imental spectra 25 ' 26,29,54 . Obviously, this kind of devia- 
tion does not simply imply the presence of a sharp mode. 
For comparison we also show in Fig. 16 (b) the spectrum 
corresponding to an unrealistic case of the normal state, 
T = 100 K, and the spin susceptibility containing the 
sharp mode (T = lOmeV, £ = 9 A). Here the onset fea- 
ture is more pronounced. 

(hi) The fine structure around 2500 cm -1 . It is sharper 
in the superconducting state than in the normal state. 
This seems not to apply to the similar structure of the 
experimental scattering rate spectra. 

(iv) The onset of the plateau of the scattering rate spec- 
tra. With increasing temperature, it shifts towards lower 
frequencies. This is consistent with the experimental 
data of Ref. 27. 

(v) The normal-state m*/m spectra — low-frequency 
asymptotics. Below ~ 300 cm , m* /m rapidly increases 
with decreasing frequency. This is an artefact of our 
computational approach which has an interesting physi- 
cal background. We calculate 02 (w) systematically using 
Eq. 26 and (K) by taking the average of the expression 
on the right hand side of Eq. 22. The value of (K) is 
very slightly (by about 2 %) larger than the one given 
by the sum rule of Eq. 27. As a consequence, o-(uj) 
contains a small unphysical singular part, which is re- 
sponsible for the behaviour of m*/m at low frequencies. 
We have checked that the problem is not due to a lim- 
ited precision of our numerical calculations. Instead, it 
seems to be caused by the neglect of vertex corrections 
in Eq. 29. The latter have been recently suggested to in- 
fluence the response functions dramatically 61 ' 62 . For the 
present values of the parameters, their role, at least as 
the integrated spectral weight is concerned, appears to be 
only minor. They can be expected, however, to become 
more pronounced for higher values of g. Indeed, we have 
found, that the discrepancy between the results of Eqs. 22 
and Eq. 27 mentioned above increases considerably with 
increasing g. Note finally that for u> > 500 cm -1 the 
spectra of m*/m are almost unaffected by the spurious 
component of a. 

(vi) High-frequency asymptotics ofm*/m. At higher 
frequencies, [m*/m}(u,T = 20K) > [m*/m](u,T = 
100 K). For the present values of the parameters, the 
two spectra converge only slowly to their common limit 
of 1. At u) = 5000 cm" 1 , the magnitude of the difference 
is still approximately 0.02. 



2. Temperature dependence of the total intraband spectral 
weight 

The changes of the total intraband spectral weight, 
Io = J °° du; (7 1 (ui), arc of particular interest, because this 
quantity is approximately proportional to the effective 
kinetic energy of the charged quasiparticles. The latter 
is defined by 

K.E. = ^2 £ k«-k , = 5Z^ c ka Cka ) ( 52 ) 

k a 

with an obvious extension to multilayer compounds. In 
the case of the nearest-neighbor tight-binding dispersion 
relation, the sum rule of Eq. 27 reduces to 

7T e 2 1 

J o = -^o x K.E./unit cell . (53) 
^ an z 

For very small values — relative to that of t — of further- 
neighbor hopping parameters, Eq. 53 can be expected to 
hold approximately. The values of the latter parameters 
characteristic of the real high-T c materials are smaller yet 
comparable with t. It is thus not immediately clear, to 
what extent Eq. 53 applies. In order to obtain an insight, 
we first concentrate on the trivial case of noninteracting 
quasiparticles (g = 0). In the remaining part of this 
subsubsection we shall discuss the TD of Io and related 
quantities, resulting from the spin-fermion model used 
above. 

Figure 17 shows the computed temperature dependen- 
cies of /1, K.E./unit cell, and Io for the dispersion rela- 
tion of Eq. 32 and of Ref. 20. It can be seen that (i) 
The chemical potential increases monotonically with in- 
creasing temperature, (ii) The kinetic energy increases 
with increasing temperature, approximately following a 
parabola, (iii) In part (a), Io is approximately equal to 
one half of — 2d?i 2 /(7re 2 )K.E., in agreement with Eq. 53, 
and the TD of Io also tracks the one of K.E. In part 
(b), the TD is more complex: Io is approximately T- 
independent for T < 200 K and it decreases at higher 
temperatures. Note that the only generic trend is that 
of the point (ii). The other aspects depend very much 
on details of the dispersion relation. Even for the simple 
t — tf dispersion — with the present values of the hopping 
parameters — both fj, and Io can either increase/decrease 
or decrease/increase depending on the location of the 
Fermi surface with respect to the Van-Hove singularity 
associated with the V-point. The first (second) variant 
occurs for the Fermi surface above (below) the singular- 
ity. Next we discuss changes of the three quantities when 
going from the normal to the superconducting state, (i) 
The chemical potential increases in both cases. This can 
be explained as follows: in the superconducting state the 
region of the BZ around the Van-Hove singularity be- 
comes partially unoccupied and the chemical potential 
has to increase in order that the particle number would 
be conserved, (ii) The kinetic energy increases, a char- 
acteristic feature of the BCS theory, (iii) In part (a) the 
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spectral weight Iq> behaves similarly as — K.E, in agree- 
ment with Eq. 53. In part (b), however, io very slightly 
increases, contrary to what one would expect based on 
Eq. 53. Again, only the trend of the point (ii) is a generic 
one. 
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FIG. 18: The same as in Fig. 17 but for quasiparticles coupled 
to the spin fluctuations. The lines correspond to the normal 
state and the solid symbols to the superconducting state. For 
each value of T the value of \i was adjusted to yield the total 
number of electrons per unit cell of 0.76. The open sym- 
bols (the symbol x and the star) correspond to an unrealistic 
case of the superconducting state, fi — /^(normal state, T = 
20 K), x (q,w) = x(normal state, T = 20K,q,w) |*(q,w) = 
X (superconducting state, T = 20K,q, u)]. 



FIG. 17: Temperature dependencies of the total intraband 
spectral weight Io, the effective kinetic energy per unit cell 
multiplied by — 2dh 2 /(ire 2 ) (labelled as — K.E.), and the chem- 
ical potential for noninteracting quasiparticles. Parts (a) and 
(b) correspond to the dispersion relations of Eq. 32 and of 
Ref. 20, respectively. For each temperature, the value of the 
chemical potential was adjusted to yield the total number of 
electrons per unit cell of 0.74 (a) and 0.84 (b). The symbols 
correspond to the BCS ansatz with the superconducting gap 
given by Eq. 33. 

Figure 18 shows the computed temperature depen- 
dencies of fi, K.E. /unit cell, and Io for quasiparticles 
coupled to the spin fluctuations. The dispersion re- 
lation of Eq. 32 and the model spin susceptibility in- 
troduced in the preceding subsection and subsubsec- 
tion have been used. The normal-state trends are sim- 
ilar to those of Fig. 17 (a): /i and K.E. increase and 
Io decreases. Both the quantities K.E., Io, and their 
changes, AI = I o (200K) - Z o (20K) and AK.E. = 
K.E.(200K) - K.E. (20 K), approximately fulfill Eq. 53. 
The changes occuring upon entering the superconduct- 
ing state (at T = 20 K) are the following: fi decreases 
by ~ 5meV, K.E. increases, and Io decreases. The 
magnitudes of the changes of K.E. and I are surpris- 
ingly similar to those of Fig. 17, A/ = io(SC,T = 
20K)-7o(N,T= 20 K) = -0.04 • 10 6 f^cm" 2 . Next 



we propose an interpretation of these effects associated 
with the superconducting transition. For the present pur- 
poses and within the context of our approach the lat- 
ter transition can be thought of as consisting of three 
stages: (i) Opening of the superconducting gap without 
relaxing the chemical potential and the spin susceptibil- 
ity, (ii) "Improving" (sharpening) of the quasiparticles 
due to the superconductivity related change of the spin 
susceptibility, in particular, due to the formation of the 
spin gap (which is, within our naive model, achieved by 
the narrowing of the low frequency component of x), still 
without relaxing /j. This is the feedback effect of the spin 
fluctuations on the quasiparticles which has been already 
considered by several authors 63,64 ' 65 , (iii) Change of the 
momentum distribution function connected with the 
change of \i which is reinforced by the requirement of 
particle-number conservation. The first stage results in 
an increase of K.E. and a related decrease of Io, see the 
open symbols in Fig. 18. These are the same effects as in 
Fig. 17 (a) and their explanation is trivial. The second 
stage is associated with a decrease of K.E. and an increase 
of Io, see the symbols + and x in Fig. 18. These effects 
are basically due to the fact that the superconducting 
state is, as regards the momentum distribution function, 
more Fermi-liquid like than the normal state 64 - 65 . If it 
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were only for the two stages, (i) and (ii), the K.E. (Iq) 
would — for the present values of the input parameters — 
decrease (increase) when going from the normal to the 
superconducting state. The third stage, however, can re- 
verse the result. This is what happens for the present val- 
ues of the input parameters. A lesson one can learn from 
the example of this subsubsection is that the changes of 
K.E. and lo associated with the superconducting transi- 
tion are determined by a delicate competition of the three 
ingredients: the opening of the superconducting gap, the 
feedback effect of the spin fluctuations on the quasipar- 
ticles, and the chemical potential shift, which may be 
large because of an inherent electron-hole asymmetry. 
So far we have considered changes occuring at a fixed 
temperature of 20 K. These are not directly accessible 
to experimental investigations and can only be obtained 
by a careful elimination of the normal state trend. The 
quantity A7 = 7 (SC,T = 20 K) - J (N,T = 100 K), 
however, can be measured. Our calculations yield a small 
and positive value of A'I of 0.015 • 10 6 Q^cm" 2 . 



V. SUMMARY AND CONCLUSIONS 

In Subsects. IVA-IVD we have analyzed the struc- 
tures of the superconducting state (SCS) conductivity 
reflecting the coupling of the charge quasiparticles to 
the resonance mode: the onset of the real part <j\ (uj) 
of the conductivity starting around the frequency ujq of 
the mode and the maximum of the function W[ui) cen- 
tered approximately at lo = ujq + Aq/H, where Ao is the 
maximum value of the superconducting gap. Our analy- 
sis provides a clear interpretation of both features. It is 
based on a division of the Brillouin into three parts: the 
nodal, the region, and the intermediate. Both <7i(w) and 
W(uu) can be expressed as the sums of contributions of 
the three parts, whose structures can be easily related to 
those of the quasiparticle spectral function A(k, to). This 
gives us a possibility to understand also the features of 
<Ji{u>) and W(ui) in terms of A(k, u). We have shown 
that the onset is due to the appearance above ujq of low 
energy excitations of the nodal region, consisting of two 
nodal quasiparticles and the mode. For the conventional 
shape of the d-gap the onset is very gradual but it be- 
comes steeper with increasing broadening of the nodes. 
The maximum of W(co) has been shown to be due to 
the appearance above ujcsb = + Aq/?i of excitations 
of the nodal region, consisting of a nodal particle, an 
antinodal particle and the mode. Due to the Brillouin 
zone averaging, in particular, due a considerable contri- 
bution of the intermediate region, the maximum of W{lo) 
appears at a frequency £Im which is slightly higher but 
considerably lower than ujq + {2Aq/K). This result of our 
analysis confirms, to some extent, the claim by Carbotte, 
Schachinger, and Basov, that the frequency of the maxi- 
mum is approximately equal to Qcsb- The width of the 
maximum, however, is determined not only by the one of 
the mode but also by the averaging. For zero width of 



the mode, the FWHM of the maximum is ca 160 cm . 

In the light of the present analysis, the interpretation 
of the onset feature in the SCS infrared spectra in terms 
of the neutron resonance remains to be very plausible. 
In fact, we are not aware of any other comparable ex- 
planation. An encouraging finding is that the values of 
the coupling constant to be used in conjunction with the 
properly normalized spin susceptibility are of the same 
order of magnitude as those of U expected to be relevant 
for the HTCS. 

In Subsect. IV E we have studied the role of the spin- 
fluctuation continuum. The computed spectra of the 
scattering rate 1/r exhibit a sharp step (onset feature) 
centered around Qm, & region of an approximately linear 
increase of 1/r, and a plateau in the mid-infrared. The 
height of the step is determined by the spectral weight 
Ir of the mode, the difference between the value of 1/r 
of the plateau and the height of the step by the spetral 
weight Ic of the continuum, and the onset frequency of 
the plateau by the width of the continuum. The exper- 
imental scattering rate spectra can be reproduced using 
values of Ic/Im of 5-10 and those of the width of several 
hundreds meV. The computed spectra of 1/r display, in 
addition, an interesting weak structure located approxi- 
mately at fi — e-p + Ao + hu>o (/z is the chemical potential 
and er the quasiparticle energy of the T-point) . The pres- 
ence of this structure in the experimental spectra would 
allow one to estimate the width of the occupied part of 
the band. 

Finally, in Subsect. IV F we have concentrated on the 
temperature dependence of the spectra. The main differ- 
ences between the SCS and the normal state (NS) exper- 
imental data, in particular, the spectral gap, the tail of 
a\(u), NS) — (Ti(w, SCS) ranging to very high frequencies, 
and the shift of the scattering rate plateau towards lower 
frequencies with increasing temperature are reasonably 
well reproduced. Some of the remaining discrepancies are 
very likely to be related to the absence of self-consistency 
and vertex corrections. Wc have further investigated the 
temperature dependence of the effective kinetic energy 
K.E. and of the intraband spectral weight lo, first for 
the BCS case and second for the present version of spin- 
fermion model. In the BCS case, K.E. increases both 
with increasing temperature in the NS and when going 
from the NS to the SCS. The spectral weight behaves sim- 
ilarly as —K.E. only for the simpler of the two (realistic) 
dispersion relations used (the "i — t'" dispersion relation) 
and for the Fermi surface above the Van-Hove singular- 
ity. For the other dispersion relation even the sign of the 
change of lo between the NS and the SCS differs from the 
one of — K.E.. The estimates of the (in-plane) kinetic en- 
ergy changes based on optical data should thus be taken 
with caution. The temperature dependencies of K.E. and 
Jo calculated for the case of the quasiparticles coupled to 
the spin fluctuations are, for the present values of the in- 
put parameters, similar to those of the BCS, t — t' case. 
The physics underlying the differences between the NS 
and the SCS, however, is more complicated. It involves 
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three different mechanisms causing changes of the occu- 
pation numbers n^: the opening of the gap [leading to 
an increase (decrease) of K.E. (Io)]- a sharpening of the 
quasiparticles in the SCS associated with the formation 
of the spin gap (resulting in an opposite trend), and a 
chemical potential shift (leading here to the same trend 
as the opening of the gap). While the mechanisms are 
certainly generic, the outcome of their competition, may 
depend on details. 
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